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ABSTRACT 

A  Nusselt-type  analysis  was  performed  for  laminar  film 
condensation  on  the  inside  of  a  rotating  truncated  cone 
with  small  half  cone  angles.   This  analysis  included  the 
interfacial  shear  between  the  vapor  and  condensate,  the 
vapor  pressure  drop,  the  thermal  resistance  in  the  condenser 
wall,  and  the  condenser  outside  cooling  mechanism.   An 
approximation  of  the  analytical  model  made  it  possible  to 
find  a  numerical  solution  for  small  half  cone  angles  greater 
than  zero. 

A  non-capillary  rotating  heat  pipe  containing  an 
evaporator,  condenser,  and  distilled  water  as  the  working 
fluid  was  tested.   It  was  rotated  at  702  and  1404  RPM ,  and 
the  heat  transfer  rates  of  the  heat  pipe  were  determined 
experimentally  for  different  saturation  temperatures  cor- 
responding to  electrical  power  inputs  ranging  from  1  kW 
to  9  kW. 

The  experimental  results  showed  that  the  non-capillary 
rotating  heat  pipe  was  an  effective  heat  transfer  device. 
The  approximate  numerical  solution  conservatively  predicted 
the  heat  transfer  rate  with  a  deviation  of  18%  at  702  RPM 
and  5%  at  1404  RPM. 
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I.   INTRODUCTION 

A.  CONVENTIONAL  HEAT  PIPF 

A  heat  pipe  is  a  closed  container  capable  of  transfer- 
ring large  amounts  of  heat  at  nearly  isothermal  conditions. 
Conventionally  it  consists   of  four  main  parts:   (1)  evap- 
orator, (2)  condenser,  (3)  working  fluid,  and  (4)  capillary 
wicking  structure  to  force  the  condensate  back  to  the 
evaporator.   The  working  fluid  is  evaporated  in  the  evap- 
orator section  of  the  heat  pipe  by  the  heat  source.   As  the 
vapor  pressure  increases,  the  vapor  flows  toward  the  con- 
denser end,  transporting  the  heat  as  the  latent  heat  of 
vaporization.   At  the  condenser  end,  the  heat  is  removed 
from  the  vapor  as  it  condenses.   The  condensate  is  then 
returned  to  the  evaporator  by  the  capillary  pumping  action 
within  the  wicking  structure. 

Studies  of  the  heat  transfer  capabilities  have  been  made 
on  the  conventional  heat  pipe  [6,  10]   with  the  conclusion 
that,  for  ordinary  fluids  and  moderate  heat  fluxes,  the  heat 
pipe  is  limited  in  its  capability  to  transfer  heat  by  the 
rate  of  condensate  that  can  be  returned  to  the  evaporator. 

B.  THE  ROTATING,  NON-CAPILLARY  HEAT  PIPF 

The  rate  of  condensate  flow  may  be  increased  by  removing 
the  wicking  structure  from  the  condenser.   If  the  condenser 


* 

Numbers  in  brackets  indicate  references  listed  in  the 
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end  is  then  internally  tapered  and  the  heat  pipe  is  rotated 
about  its  longitudinal  axis,  a  centrifugal  force  field  is 
set  up  which  will  overcome  both  the  viscous  losses  in  the 
system  and  the  gravitational  field. 

Ballback's  study  [10]  of  the  rotating  heat  pipe  indi- 
cated that  its  heat  transfer  capability  is  limited  by  the 
amount  of  heat  which  can  be  removed  from  the  condenser.   In 
his  study,  Ballback  performed  a  Nusselt-type  analysis  for 
the  condensation  on  the  inside  of  a  truncated  cone.   In 
order  to  obtain  a  closed  form  approximate  solution  for  the 
film  thickness,  6,  he  made  the  following  assumptions:   (1) 
-5—  «  tan  (J)  (large  cone  angles)  ,  (2)  no  shear  stresses  between 
the  vapor  and  condensate,  and  (3)  no  thermal  resistance  in 
the  condenser  wall.   Daley  [11]  was  able  to  improve  on  the 
theory  introduced  by  Ballback  by  taking  into  consideration 
the  thermal  resistance  in  the  condenser  wall  and  the  outside 
cooling  mechanism.   Furthermore,  Daley  neglected  the  assump- 
tion that  -t—  «  tan  (J).   Using  the  velocity  profile  for  the 
condensate  developed  by  Ballback,  he  arrived  at  a  non-linear 
second  order  differential  equation  for  the  film  thickness, 
6 .   Daley  numerically  integrated  this  expression  using  a 
Runge-Kutta-Gill  numerical  integration  scheme  with  an 
IBM  360  MOD  67  digital  computer.   To  start  the  integration, 

he  used  the  conditions  that  6  =  6i  and  3—  =  tand)  at  X  =  0 . 

dx       Y 

It  was  postulated  that  6i  was  a  function  of  the  over-fall 
condition  at  the  step  located  at  the  condenser  exit.   Daley 
arrived  at  an  expression  for  the  minimum  film  thickness,  6m, 


at  the  condenser  exit,  using  the  approximation  that  the  flow 
over  the  step  approximated  the  free  fall  condition  for  open 
channel  flow.   Having  found  the  film  thickness  from  the 
second  order  differential  equation  at  X  =  Lc,  the  condenser 
exit,  Daley  compared  this  value  with  6m  at  Lc.   When  the 
difference  was  less  than  .004,  a  solution  was  found.   Using 
this  technique,  Daley  was  able  to  find  the  heat  flow  out  of 
the  condenser  as  a  function  of  RPM  for  0°  half  cone  angle 
only.   As  the  cone  angle  was  increased,  however,  instability 
of  the  highly  non-linear  second  order  differential  equation 
for  the  condensate  film  thickness  caused  the  equation  to  be 
extremely  sensitive  to  the  initial  value  of  the  film  thick- 
ness.  Consequently,  no  solutions  were  found  for  the  conden- 
sate film  thickness  for  half  cone  angles  greater  than  zero. 

C.   THESIS  OBJECTIVES 

The  objectives  of  this  study  were  to:   (1)  continue  the 
theoretical  study  of  the  heat  transfer  capability  of  the 
rotating  heat  pipe  and  obtain  solutions  for  $  >  0° ,  (2) 
obtain  experimental  data  from  the  heat  pipe  designed  and 
assembled  during  Daley's  study,  and  (3)  compare  the  analyti- 
cal and  experimental  results. 
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II.   THEORETICAL  PROGRAM 

In  order  to  improve  on  the  analysis  performed  by 
Ballback  and  Daley  the  effects  of  vapor  pressure  drop  and 
interfacial  shear  stress  between  the  liquid  and  the  vapor 
were  included.   Furthermore,  the  thermal  resistance  in  the 
condenser  wall  and  the  outside  cooling  mechanism  were  con- 
sidered.  Using  the  coordinate  system  in  Figure  1,  for 
small  half  cone  angles,  momentum,  continuity,  and  energy 
equations  were  developed  for  an  infinitesimal  fluid  element 
In  addition,  assuming  one  dimensional  turbulent  flow  in  the 
vapor,  momentum  and  continuity  equations  were  developed  for 
the  isothermal  vapor  flow. 


Figure  1.   Coordinate  system  and  geometry 
for  condensate  film  analysis. 
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A.  CONDENSATE  MOMENTUM  EQUATION  (X~Direction) 

By  assuming  that  the  momentum  changes  in  the  condensate 
are  small  a  static  force  balance  may  be  taken  on  an  infini- 
tesimal liquid  element  in  the  x-direction  as  shown  in 
Figure  2 . 

0  =  57"  §f+  Of"2  rsi**  U> 

where 

t  =  shear  stress 

p  =  pressure 

x  =  co-ordinate  measuring  distance  along  surface 

y  =  co-ordinate  measuring  distance  normal  to  surface 

r  =  radius 

<f)  =  half  cone  angle 

pf  =  density  of  liquid 

to  =  angular  velocity 

B.  CONDENSATE  MOMENTUM  EQUATION  (Y-Direction) 
Similarly,  using  Figure  2,  a  force  balance  in  the 

Y-direction  may  be  taken 

0  =  -  |£  -  p.w2  r  cost})  (2) 

dy     r 

Making  use  of  Figure  2,  the  radius  at  any  point  in  the  fluid 

condensate  may  be  expressed  as : 

r(x,y)  =  R  +  x  sin  <t>  -  y  cos  $  (3) 

where 

R  =  minimum  wall  radius  in  the  condenser  section 
o 
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Figure  2.   Cross-section  of  infinitesimal  liquid 
element. 
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Substituting  equation  (3)  for  r(x,y)  into  (2),  and  inte- 
grating (2)  from  p  to  p   and  y  to  6  yields: 


p(x,y)  =  Pv(x)  +  pfco  cos  <J){(Ro+ x  sin(f))  (6-y) 


5-  cos<J)(62-y2)} 


(4) 


Differentiating   equation    (4)    with   respect   to   x   gives 


3^-  =   g~  +    Pfw   cos  (}>{(R0+  x  sin  4>-6cos  <J>)    ^ 


+    (6-y)    sin    <j>} 


(5) 


Substituting  equation  (5)  into  equation  (1)  and  integrating 
from  y  to  6  and  x  to  -x   gives 


x  = 


dx 


(6-y) 


p. 00 


'cos  <J){(R  +x  sin  (J)- 6  cos  (J))  ^~(6-y)} 


+  p^w  sin(j){R  +xsin  (J)  -  6  cos  (J)}(6-y) 


=  y 


9u(x,y) 


f  9y 


(6) 


where 

uf     =  viscosity  of  the  liquid 
u(x,y)  =  velocity  of  the  liquid 

Integrating  equation  (6)  from  0  to  y  and  from  0  to  u  gives 
an  expression  for  the  velocity  of  the  fluid: 


u(x'y)    =    "    y^  dx~    (6y    "      2  )    "    y7" 


p.or 


f    L 


R  +  x  sinc{>-6coscj) 


sine}) -cos  (J) 


21 


i5y-  J 


(7) 


14 


This  expression  for  u(x,y)  agrees  with  Ballback's  except 
for  the  shear  stress  term. 

C.   CONTINUITY  EQUATION 

It  is  apparent  from  Figure  3  that  due  to  the  conserva- 
tion of  mass,  with  condensation 


dw.e    dw 
f  _  v 

dx    dx 


(10) 


wf  =  wv         .  (10a) 

where 

wf  =  mass  flow  rate  of  fluid 
w   =  mass  flow  rate  of  vapor 

The  mass  flow  rate  of  the  liquid  may  be  written  as 

6 
wf  =  /   Pfu  2^r  dy  (11) 

o 

The  integrating  of  equation  (11)  using  equations  (3)  and 
(7)  yields 

wf  =  2irpf  {-   ~  ^  [(RQ+  x  sin  (j>)    —.  -    cos  <J>  ^  6    ]} 

T                                     fi2  .3 

+    27Tpf{-   — -[  (Rq+  x  sine}))    -^ cos  $-3-]} 

2    2 

27Tpf      W  -jr 

+   (R  +xsind>-6  costf>)(sin  $  -  cos  §  -5—) 

u        o  dx 

[(Rq  +  xsin^)y  -  cos  d-  2J  <$4  ]  (12) 

The  mass  1  ate  of  flow  of  the  vapor  may  also  be  written 
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w     +-^ <3x 

v      cjx 

,+  ^ax)  *  +  •£&>* 

dx  v 


dx 

dp, 
V    dx" 


dx- 


7 7 7"" 

dx   _ 


dw.p 

wf+  s  dx 


V- 


Fiaure    3.      Cross    section   of    infinitesimal 
Figure  liquid  and  vapor  elements. 
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wv   =    PvttR2v  (13) 

where 

R  =   R     +  x  sin  d)  -  6cos  A 

O  T  Y 


D.  ENERGY  EQUATION  IN  CONDENSATE 

From  Figure  2,  using  an  energy  balance,  the  rate  of 
condensation  may  be  written 

,n  dwf         2tt(R  +x  sin  <J>)  (T    -T    ) 

dx-  -  hfg  -^ XTXTT (14) 

k.        k  h 

f  w 

where 

T   =  saturation  temperature  of  vapor 

Tm  =  inlet  coolant  temperature 

t   =  thickness  of  condenser  wall 

6   =  film  thickness  of  condensate 

k_  =  thermal  conductivity  of  condensate 

k   =  thermal  conductivity  of  condenser  wall 
w  2 

h   =  outside  condenser  wall  heat  transfer 
coefficient  , 

E.  VAPOR  MOMENTUM  EQUATION 

Using  Figure  3,  and  neglecting  the  momentum  of  the 
condensing  vapor,  a  force  balance  may  be  taken  in  the 
x-direction. 

dPV 
PvA  -  (pv+dpv)  (A-dA)  -  (pv+-T-)dA  -  TvdAw 

=  (w  +dw  ) (v+dv)  -  w  v  (15) 

where 
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Pv  -   pressure  of  vapor 

v   =  velocity  of  vapor 

A  =  cross  sectional  area  of  vapor  element 

Aw  =  surface  area  over  which  the  shear 
stress  acts 

Neglecting  higher  order  terms,  and  using  the  fact  that 

dA^  =  2irRdx  equation  (15)  reduces  to 

TTK 

1     2 
Using  the  expression  tv  =  -  pvv  f,  where  f  =  friction  factor, 

and  substituting  into  equation  (16)  yields 

dPv         dv    v   d  ,     7    ,     L  pvv2f 

sr  =  ~  pvv  &  -  ^2  ^  (vTR  v)  +  -5-         (17) 

Expanding  the  derivative  in  equation  (17)  yields 

*S  .  _  2p  v  dv  .  fAr  dR  +  ^  (17a) 

dx        v   dx     R   dx      R 

The  first  two  terms  in  equation  (17)  represent  the  momentum 
effects  of  the  vapor  while'  the  last  term  represents  the 
usual  friction  effect.   For  the  low  condensate  velocities 
which  were  expected  during  the  experimental  runs ,  according 
to  Bergelin  et  al.  [2]  the  friction  factor  is  approximately 
equal  to  the  friction  factor  used  in  smooth  tubes  (that  is, 
the  effect  of  the  liquid  film  on  the  condenser  surface  is 
negligible).   According  to  Rohsenow  and  Choi  [3],  the 
friction  factor  for  turbulent  flow  in  a  smooth  tube  for  vapor 
Reynolds  numbers  between  2000  and  50,000  is  expressed  as: 
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0791,     ,    .0791, 
7 (Re) *  = 


m 


(18) 


F.   DETERMINATION  OF  FILM  THICKNESS  AND  HEAT  TRANSFER  RATE 

Equations  (10) ,  (12) ,  (13) ,  (14) ,  and  (17a)  can  be 

further  reduced  to  two  simultaneous  ordinary  first  order 

differential  equations.   Using  equation  (14)  with  the  use 

dWf 
of  equations  (10)  and  (13)  for yields: 


dv    ^<V  X  Sin  ®   (Ts  -  T^    2v  r 
d>: 


kf   kw 


R2h 


fc 


sin©  -  cos 


d6 

dx 


(19) 


At  the  same  time,  using  equation  (12)  in  conjunction  with 
(10a)  and  (13)  we  obtain 


Pv7TR2V    =     2TTPf|-     J-    _ 


UK   +  x  sirnji) 


3        C°5  y  24 


5     ,4, 


62  63 

+    2iTpf  |--p  [  (Rq+  x  sin<f>)    -J cos  <J>  -y] 


2    2 

27Tpf     0)  d6 

+    (R  +  x  sin<J>-  6cos(f))  (sin(j>-  cos    — ) 


dx 


.  [  (R  +  x  sin4>)    — 


a     5     x4    , 
cos  <J>  2J  6       J 


(20) 


If    equation    (17a)    is   now   substituted   into   equation    (20) ,    an 
equation    for   —  is    obtained: 


d6_ 
dx 


dx 


tan 


Pvv2f 


R         2 


+    Pvv 


vfR' 

~2p; 


-2p 


dv 

1   dx 


p.  cos 


2         2pvv 
pfW    R+-F- 


where 


(R  +  x  sin  <P)    -3 cos  <J>  24-  6 


p      =    (R  +  xsin<J>)    -5 cos  <J> 


63 


(21) 
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Equations  (19)  and  (21)  are  two  simultaneous  first 
order  ordinary  differential  equations  with  unknowns  6 (x) 
and  v(x).   A  Runge-Kutta  numerical  integration  was  used  to 
solve  them  with  the  IBM  digital  computer.   It  was  hoped 
that  by  using  first  order  differential  equations  rather 
than  a  second  order,  as  Daley  did,  that  the  sensitivity  to 
initial  6  would  be  reduced.   Figure  4  is  a  flow  chart  of 
the  program  used  to  solve  the  differential  equations. 

The  integration  was  started  by  assuming  initial  con- 
ditions, 6=6.  and  v=0  at  x=0.   With  these  initial  conditions 
-3—  =  tan  4)  at  x=0.   Figure  5  shows  a  schematic  representing 
the  calculated  6 (x)  for  different  initial  values  of  6. 
Note  that  if  6.    is  too  small  the  calculated  5  becomes 
negative  before  the  condenser  overfall  (physically  impossible 
solution) .   If  6 .  is  too  large  the  calculated  6  grows  at  a 
constant  slope  equal  to  tan<f>  (physically  unrealistic).   With 
the  proper  choice  of  6 .  the  calculated  6  reaches  a  maximum 
along  the  length  of  the  condenser  and  then  tapers  off  to  a 
value  close  to  zero  at  the  overfall.   The  exact  value  of  6 
at  the  overfall  is  not  known.   It  was  observed  however  that 
the  calculated  value  of  6  at  the  overfall  was  approximately 

.1  of  6.,  and  the  heat  transfer  rate  was  insensitive  to  an 
l 

exact  value  at  the  overfall.   This  result  agrees  with  Nimmo 
and  Leppert  [8]  who  experimentally  found  for  laminar  film 
condensation  on  a  finite  horizontal  surface  that  the  heat 
transfer  rates  were  practically  independent  of  the  end 
conditions . 
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Q   START   J 


)Z 

INPUT  DATA 

V  pv,  Pf,  uf,vv,kf 

hfg'  V  h'  Lc'  V 
<f>,  T  ,  t,  RPM 


V 


INITIAL  CONDITIONS 
6=6,  ^ 


v=0 


x=0 


Ai 


CALCULATION  OF  Q 

Using  the  value  of  6  found  in 
Subroutine  RK1  equation  (14)  is 
used  to  find  the  heat  transfer 
rate  at  each  interval  of  x.   Q 
is  obtained  by  using  the  trap- 
ezoidal rule. 


*r 


SUBROUTINE  RKl 

Numerically  integrates 
equations  (19)  and  (21) 
using  1000  intervals 
of  x  from  x=0  to  x=L 


PRINT  6,    v,    Q 


7V 


t 


SUBROUTINE  DFRS 

Contains  explicit  func- 
tions of  derivatives, 

dv    ,  d6 
-s—  and  3—  . 
dx      dx 


Figure  4.   Flow  chart  of  computer  program. 
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tan 


2  Tpl   p?~ 


+  p  v 


rvfR' 

l2"p7 


2pi  dT 


cos  <f>  p^oj  R  + 


2p  v' 

1  f  R 

dv  =  2(R0+xsin^)(Ts-Tj 


2  v  sine 


(22) 


(23) 


Equation  (22)  may  be  further  simplified  by  an  order  of  mag- 
nitude analysis.   Neglecting  the  second  terms  of  p,  and  p 


and  using  the  expression,  R  =  R  +  x  sin<j>   with  the 

av   yfR2    p-,   p9 

previously  observed  result  that  p,  ~<<  -± and  —  «  -£-, 

cl    dx    2pf       R     2 

equation  (22)  may  be  rearranged  in  terms  of  a  cubic  equation 
for  6 (x) : 


R  sin< 


2 

pfoj  R 


2pvV 


p  v2f 
^v 


2i 


yfR' 


pvV  2p 


=  0 


(24) 


Equations  (23)  and  (24)  were  solved  simultaneously  for 
v(x)  and  6  (x)  .   Equation  (-23)  was  numerically  integrated 
using  1000  steps  of  integration  over  the  cooling  length  of 
the  condenser.   The  integration  was  started  with  initial 
conditions  6.  =  0  and  v  =  0  at  x  =  0.   Having  found  v(x)  at  a 
particular  step,  equation  (24)  was  solved  by  a  subroutine. 
This  value  of  6  was  used  to  find  v(x)  for  the  next  step  of 
integration.   The  results  of  approximation  1  are  recorded 
in  Table  I  for  half  cone  angles  of  .1°,  .2°,  and  .3°. 

Using  approximation  2  on  equation  (21),  equation  (24) 

was  further  simplified  and  rearranged  to  yield  the  following 

expression  for  6 (x) : 
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TABLE  I 

THE  HEAT  FLOW,  Q  (Btu/hr) ,  FOR 

DIFFERENT  HALF  CONE  ANGLES 


NUMERICAL  APPROXIMATE 

SOLUTION  SOLUTION  1 

Eqs.  (19)  &  (21)    Eqs.  (23)  &  (24) 


APPROXIMATE 
SOLUTION  2 
Eqs.  (23)  &  (25) 


0° 

8,505 

1° 

10,082 

2° 

11,173 

3° 

11,904 

9,836 
11,033 
11,769 


10,290 
11,332 
11,998 


<5(x) 


3pvyfv 


-   2  2  . 
2pf  a)  sin 


1/3 


(25) 


Using  equations  (23)  and  (25)  a  numerical  solution  was  found 
for  v(x)  and  6 (x)  in  the  same  manner  as  the  solution  for 
approximation  1.   The  results  are  tabulated   in  Table  I  for 
half  cone  angles  of  .1°,  .2°,  and  .3°. 

Table  I  shows  that  approximation  1  is  lower  than  2. 
This  is  to  be  expected  since  approximation  1  includes  both 
interfacial  shear  between  the'  vapor  and  condensate  and  the 
pressure  drop  in  the  vapor.   At  the  same  time  approximation 
1  is  less  than  the  more  accurate  solution  of  equations  (19) 
and  (21) ,  thus  giving  a  conservative  result  with  a  deviation 


*In  finding  the  solution  for  .3°,  10,000  steps  were 
used  in  the  integration.   This  was  necessary  because  as  the 
half  cone  angle  became  larger  than  .2°,  the _ numerical  inte- 
gration became  very  sensitive  to  the  step  size. 
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of  only  2.5%.   Consequently,  approximation  1  v/as  used  to 
arrive  at  a  theoretical  solution  for  6 (x)  and  Q  for  the 
one  degree  half  cone  angle  condenser  used  in  the  experi- 
mental, program.   The  results  are  plotted  in  Figures  9 
and  10. 
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II.   EXPERIMENTAL  PROGRAM 

A.   DESCRIPTION  OF  EQUIPMENT 

In  order  to  facilitate  the  description  of  the  entire 
experimental  system,  the  system  is  sub-grouped  into  four 
parts:   (1)  evaporator,  (2)  condenser,  (3)  instrumentation, 
and  (4)  auxiliary  equipment.   Included  in  the  description 
are  drawings  which  are  intended  to  aid  the  reader  in 
visualizing  the  different  components  of  the  system. 
Figure  6  is  a  schematic  representation  of  the  entire  system 
including  the  auxiliary  equipment.   Figure  7  is  a  cross- 
sectional  drawing  of  the  assembled  heat  pipe.   Figure  8 
is  a  photograph  of  the  heat  pipe  and  associated  equipment. 

1.   Evaporator 

The  evaporator,  modeled  to  approximate  the  rotating 
boiler  tested  at  Lewis  Research  Center  [17] ,  is  made  of  304 
stainless  steel.   It  has  an  inside  diameter  of  3.125  inches, 
an  outside  diameter  of  3.635  inches,  an  exit  diameter  of 
2.625  inches,  and  an  overall  length  of  5.905  inches. 

The  actual  heated  length  is  3-7/8  inches  around  which 
an  11-gauge  Chromel-A  heater  wire  in  a  3/16  inch  inconel 
sheath  is  helically  wound.   The  procedure  of  bonding  the 
heater  wire  to  the  evaporator  is  an  exacting  process  in 
order  to  ensure  a  uniform  heat  flux  axially  along  the  heater 
surface  and  to  guard  against  electrical  shorting  between 
the  heater  wire  and  other  components.   The  procedure  was  as 
follows : 
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a.  A  sheet  of  gold-nickel  foil,  .003  inch  thick, 
was  tack  welded  to  the  outside  evaporator  surface. 

b.  The  heater  wire  was  helically  wound  over  the 
gold-nickel  sheet  with  two  wraps  of  .020  inch  diameter  gold- 
nickel  wire  wound  intermittently  between  the  heater  coils. 

c.  The  assembly  was  then  brazed  in  a  hydrogen 
atmospherically-controlled  oven . 

d.  Sauereisen  cement  was  baked  around  the  outside 
of  the  heater  wire  using  the  heater  wire  as  the  heat  source, 

e.  Asbestos  insulation  was  wrapped  around  the 
heater  wire  and  tied  off  with  glass  tape. 

f.  Fish  paper  was  placed  into  the  heat  sink  on  the 
condenser  end  of  the  evaporator  to  insulate  the  heater 
coils  from  the  evaporator. 

Power  was  furnished  to  the  heater  coils  by  a  DC  motor- 
generator  capable  of  delivering  150  amperes  at  220  volts. 
The  electrical  power  was  passed  by  a  brush  assembly  con- 
sisting of  eight  spring-loaded  graphite  brushes  through  two 
bronze  collector  rings  to  the  helically  wound  heater  wire. 
The  rings  were  supported  by  twelve  1/4  inch  stud  bolts. 
Phenolic  spacers  were  used  to  align  the  collector  rings  and 
to  insulate  the  bolts  from  the  rings. 

In  order  to  reduce  the  axial  heat  losses,  two  radial 
1/16  inch  wide  grooves  were  machined  in  the  outside  surface 
to  within  1/16  inch  of  the  inside  boiler  wall  at  either  end 
of  the  heating  element. 
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The  outside  face  of  the  evaporator  is  closed  by  the 
insertion  of  two  pyrex  glass  view  ports,  3-1/2  inches  in 
diameter.   The  inner  view  port  rests  on  a  neoprene  O-ring 
seal  and  is  separated  from  the  outer  view  port  by  a 
3/32  inch  compressed  fiber  gasket.   The  view  ports  are  held 
in  place  by  a  stainless  steel  end  flange  secured  in  place 
by  twelve  1/4  inch  stud  bolts.   The  original  design  of  two 
solid  glass  view  ports  was  changed  by  fusing  the  inner  view 
port  to  a  glass  tube  which  passes  through  the  outer  view 
port.   This  tube  allows  a  vacuum  to  be  drawn  from  the  inside 
of  the  heat  pipe  during  the  filling  procedure,  and  is  easily 
sealed. 

The  inside  diameter  of  the  evaporator  is  stepped  1/2 
inch  from  the  exit  diameter  to  allow  for  the  formation  of  a 
liquid  annulus  in  the  evaporator. 

2 .   Condenser 

The  condenser  section  is  machined  of  stainless  steel 
in  the  shape  of  a  truncated  cone.   The  half  cone  angle  is 
1°.   The  base  of  the  condenser  mates  to  the  evaporator  exit 
and  has  an  inside  diameter  of  1.809  inches  and  tapers  over 
a  10  inch  length  to  a  diameter  of  1.460  inches  at  the 
truncated  end.   At  this  point,  the  condenser  extends  an 
additional  1-1/2  inches  to  house  a  cylindrical  end  plug. 

A  step  is  machined  at  the  evaporator-condenser  inter- 
section where  the  condensate  flows  from  the  condenser  to  the 
recessed  inside  diameter  of  the  evaporator.   The  first  inch 
of  the  condenser  surface  near  the  evaporator  is  assumed  to 
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be  adiabatic  and  9  inches  of  the  condenser's  overall  length 
is  considered  as  the  cooling  surface.   The  wall  thickness 
of  the  condenser  is  1/16  inch  in  order  to  reduce  its  thermal 
resistance. 

The  condenser  section  has  two  flanges.   The  large  flange 
mates  to  the  bearing  seat  flange  on  the  evaporator  and  is 
held  in  place  by  the  1/4  inch  stud  bolts.   The  other  flange 
mates  to  a  flange  on  the  end  plug  held  together  by  eight 
1/8  inch  stud  bolts.   The  end  plug  has  an  extension  on  its 
open  end  on  which  are  mounted  four  phenolic  junction  boards. 

3.   Instrumentation 

A  schematic  drawing  of  the  instrumentation  set-up 
is  shown  in  Figure  6. 

a.   Temperature  Measurements. 

Temperature  measurements  are  obtained  from  nine 
grounded  copper-constantan  thermocouples  and  two  quartz 
thermometers.   The  thermometers  are  used  to  measure  a 
temperature  profile  in  the  evaporator  wall,  the  liquid 
annulus  temperature  in  the  evaporator,  the  vapor  saturation 
temperature  in  the  evaporator,  the  outside  wall  temperature 
of  the  condenser,  and.  the  temperature  at  the  location  of  the 
pressure  transducer. 

In  order  to  measure  the  temperature  gradient  in  the 
evaporator  wall,  four  thermocouples  with  a  1/16  inch  dia- 
meter stainless  steel  sheath  are  utilized.   Fach  thermo- 
couple is  placed  within  a  1/16  inch  well  drilled  3-1/2 
inches  deep  into  the  evaporator  wall.   The  thermocouple 
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junction  is  coated  with  epoxy  to  give  it  enough  strength  so 
that  it  will  not  bend  when  placed  against  the  end  of  the 
well.   The  epoxy  is  sanded  such  that  the  junction  makes 
contact  with  the  wall. 

One  thermocouple  is  located  1/16  inch  radially  from  the 
inside  of  the  evaporator,  another  is  located  1/16  inch 
radially  in  from  the  outside  evaporator  surface,  and  the 
remaining  two  are  located  opposite  one  another  in  the  middle 
of  the  1/4  inch  wall  at  the  same  radial  position  to  check 
on  the  uniformity  of  the  temperature  field. 

The  leads  from  all  of  the  thermocouples  in  the  evapor- 
ator wall  are  brought  through  slots  in  the  evaporator  and 
forward  condenser  flanges,  down  the  side  of  the  condenser, 
and  through  the  rear  condenser  flange  to  phenolic  junction 
boards.   From  the  boards,  the  leads  go  inside  the  shaft  and 
down  to  the  slip  ring  unit  which  is  attached  to  the  main 
assembly  by  a  coupling  mounted  at  the  end  of  the  shaft. 

A  thermocouple  is  suspended  in  the  liquid  annulus  to 
measure  the  liquid  bath  temperature.   The  remaining  thermo- 
couple in  the  evaporator  is  suspended  in  the  vapor  space  to 
measure  the  saturation  temperature. 

Since  these  two  thermocouples  protrude  into  the  inside 
of  the  heat  pipe  where  a  vacuum  is  held,  precautions  were 
taken  to  ensure  that  a  seal  was  maintained.   The  ends  of 
the  two  thermocouples  were  enclosed  in  stainless  steel 
flanges  and  the  leads  were  soldered  at  the  slots  on  the 
condenser  forward  flange  to  keep  air  from  passing  through 
the  slots  into  the  evaporator. 
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In  order  to  measure  the  outside  condenser  wall  temper- 
ature, two  thermocouples  were  soldered  to  grooves  on  the 
outside  of  the  condenser  at  locations  of  1/3  and  2/3  of  its 
length  respectively. 

The  remaining  thermocouple  was  used  to  measure  the 
temperature  at  the  pressure  transducer.   The  junction  is 
soldered  to  the  inside  of  the  pressure  transducer  arm  1/32 
inch  from  the  surface. 

The  temperature  of  the  coolant  is  measured  before  and 
after  spraying  by  quartz  thermometers  mounted  in  mixing 
tubes  located  above  and  below  the  condenser  spray-  cooling 
cylinder. 

b.  Pressure  Measurements 

The  saturation  pressure  of  the  vapor  at  the 
small  end  of  the  condenser  is  measured  by  a  semi-conductor 
pressure  transducer  mounted  in  the  vapor  space.   The  trans- 
ducer is  temperature  compensated  and  is  threaded  on  a  1/4 
inch  diameter  arm  that  extends  from  the  shaft  through  the 
end  plug  and  about  three  inches  into  the  vapor  space.   The 
leads  of  the  probe  pass  through  the  transducer  arm  up 
through  the  shaft  to  the  junction  boards.   From  there,  the 
leads  pass  back  through  the  shaft  to  the  slip  ring  unit. 

c.  Coolant  Flow  Rate  Measurements 

The  coolant  flow  rate  is  measured  by  passing 
the  coolant  through  a  flow  rotameter  before  the  coolant 
enters  the  spray  nozzles. 
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d.  Rotational  Speed  Measurements 

The  rotational  speed  of  the  heat  pipe  is 
measured  by  two  meters.   A  Beckman  Eput  Meter  measures  the 
RPM  using  a  magnetic  pickup,  and  an  electronic  strobe  scope 
is  used  as  a  check. 

e.  Electrical  Povzer  Measurements 

Power  to  the  heater  coils  is  measured  by  an 
ammeter  and  voltmeter  connected  in  series  and  parallel 
respectively  with  the  power  lines  returning  from  the  heater, 
The  power  level  is  controlled  by  a  field-rheostat  placed  in 
parallel  with  the  DC  motor-generator.   A  150  ampere  circuit 
breaker  is  included  in  the  circuit. 

f.  Slip  Ring  Unit 

The  slip  ring  unit  allows  the  rotating  thermo- 
couple leads  and  pressure  transducer  leads  to  pass  electric 
current  to  stationary  leads.   The  slip  ring  unit  used  in 
the  experiment  utilizes  liquid  mercury  as  a  viscous  contact 
medium  between  the  rotating  and  stationary  leads.   The 
frictional  heating  noise  is  less  than  10  microvolts  at 
speeds  up  to  4000  RPM.   The  slip  ring  unit  has  22  terminals 
to  accommodate  the  nine  thermocouples  and  pressure  trans- 
ducer.  The  output  from  the  slip  rings  is  fed  to  a  Hewlett- 
Packard  2010C  Data  Acquisition  System  containing  an  inte- 
grating digital  voltmeter,  guarded  data  amplifier,  and 
guarded  cross-bar  scanner.   This  system  provides  printed  or 
visual  output  with  an  accuracy  of  ±  .5  microvolts. 
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4 .   Auxiliary  Equipment 

The  auxiliary  equipment  includes:   (a)  drive 
assembly,  (b)  test  stand,  (c)  safety  shields,  and  (d)  spray 
assembly . 

a.   Drive  Assembly 

The  drive  assembly  consists  of  a  drive  motor, 
end  plug,  shaft,  pulley  and  bearing.   The  drive  motor  drives 
the  shaft  by  means  of  a  2.65  inch  diameter  pulley  and 
V-belt.   The  driving  motor  is  a  2  horsepower,  3-phase 
variable  speed  motor  capable  of  a  maximum  speed  of  4000  RPM. 
The  speed  is  controlled  by  a  pushbutton  electric  remote 
control  unit  and  the  motor  can  be  shut  down  instantaneously 
by  a  magnetic  disc  brake. 

The  end  plug  has  a  3/4  inch  threaded  shaftway  machined 
into  it.   Between  the  end  plug  and  the  shaftway  is  a  .27 
inch  hole  to  allow  for  the  insertion  of  the  pressure 
transducer.   The  stainless  steel  shaft  is  hollow  to  allow 
for  the  passage  of  the  thermocouple  and( the  transducer  leads 
to  the  slip  ring  unit.   The  shaft  has  3/4  inch  of  its 
outside  diameter  threaded  for  screwing  into  the  shaftway, 
and  is  threaded  for  3/8  inch  along  its  inside  diameter  on 
both  ends.   The  condenser  end  mates  with  the  transducer  arm 
described  earlier  and  the  opposite  end  is  joined  to  a 
coupling  connected  to  the  slip  ring  unit. 

There  are  two  main  bearings  supporting  the  heat  pipe. 
The  large,  5.96  inch  diameter,  floating,  single  bearing 
seats  on  the  evaporator  flange  which  is  supported  by  an 


36 


aluminum  pillow  block.   Copper  cooling  lines  through  which 
tap  water  flows  have  been  placed  in  the  pillow  block  to 
help  cool  the  bearing.   There  is  an  oil  input  line  at  the 
top  of  the  pillow  block  which  allows  for  lubrication  of  the 
bearing.   The  small  double  row,  angular  contact  bearing  on 
the  condenser  end  is  mounted  on  a  seat  machined  on  the  shaft 
and  is  self-lubricated. 

b.  Test  Stand 

The  test  stand  was  designed  and  built  so  that 
the  heat  pipe  could  be  rigidly  supported  with  its  longi- 
tudinal axis  ranging  in  position  from  0  to  90  degrees  from 
the  horizontal.   Both  the  heat  pipe  and  variable  drive 
motor  are  mounted  on  steel  bed-plates  connected  by  a  4  inch 
iron  pipe.   The  iron  pipe  is  held  in  place  by  three  2-inch 
thick  clamps  which  are  supported  two  feet  off  the  ground 
by  a  1/4  inch  steel  plate.   As  the  heat  pipe  is  rotated  in 
any  position,  the  drive  motor  rotates  with  it  and  the  entire 
assembly  is  then  securely  held  in  place.   The  motor  has  been 
vibration  mounted  in  order  to  keep  the  vibration  to  a 
minimum. 

c.  Spray  Cooling  Assembly 

Cooling  of  the  condenser  section  is  done  by 
spraying  tap  water  on  the  outside  of  the  condenser  during 
rotation.   The  condenser  section  is  enclosed  by  a  9-inch 
long  13-inch  diameter  stainless  steel  cylinder,  the  con- 
denser box.   The  ends  of  the  condenser  box  have  Garlock 
rubber  seals  separating  them  from  the  condenser  section. 
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The  condenser  box  is  slid  into  place  over  the  condenser 
section  through  its  open  ends  and  is  then  bolted  in  place 
to  the  support  plate.   An  8-inch  x  6-inch  window  has  been 
cut  in  the  condenser  box  in  order  to  be  able  to  make 
adjustments  or  repairs  without  removing  the  condenser  box. 
During  operation,  a  steel  plate  is  bolted  over  the  window. 
The  condenser  box  has  eight  spray  nozzles  mounted  at  the 
same  axial  position.   The  nozzles  may  be  varied  to  give  a 
range  of  droplet  sizes  from  300  to  600  microns. 

The  coolant  flows  from  copper  feed  lines  to  a  mixing 
tube  located  at  the  top  of  the  condenser  box.   A  quartz 
thermometer  records  the  temperature  of  the  coolant  in  the 
tube.   From  the  tube,  the  coolant  flows  through  plastic 
tubing  to  each  of  the  spray  nozzles.   The  coolant  drains 
from  the  bottom  of  the  condenser  box  to  another  mixing  tube 
where  the  temperature  is  recorded  by  the  other  quartz 
thermometer.   From  there  the  coolant  flows  to  an  external 
drain . 

d.   Safety  Shields' 

During  operation,  the  entire  heat  pipe  assembly 
is  surrounded  by  1/8  inch  thick  stainless  steel  shielding 
bolted  to  the  support  plate  to  ensure  safe  operation. 

The  evaporator  shield,  shaped  like  an  inverted  U,  is 
bolted  to  a  1.0  inch  thick  plexi-glass  window  on  its  open 
end.   This  window  serves  as  a  safety  shield  and  also  allows 
the  internal  pipe  operation  to  be  observed.   The  other  end 
mates  with  the  condenser  box. 


B.   EXPERIMENTAL  PROCEDURES 

The  experimental  procedures  consist  of:   (1)  Filling 
Procedure  and  (2)  Operational  Procedure. 

1.   Filling  Procedure 

In  order  to  obtain  film  condensation  in  the  conden- 
ser, the  inside  surface  of  the  condenser  must  be  prepared 
so  that  proper  "wetting"  will  take  place  during  operation. 
The  following  procedure  is  used: 

a.  Scrub  the  condenser  surface  with  190 °F  Alconox 
detergent  solution. 

b.  Flush  with  distilled  water. 

c.  Spray  with  alcohol. 

d.  Rinse  with  distilled  water. 

e.  Spray  with  acetone. 

f .  Rinse  with  distilled  water. 

Having  completed  the  above  procedure,  the  heat  pipe  is 
ready  to  be  filled.   In  this  study,  distilled  water  was  used 
as  the  working  fluid. 

In  order  to  get  a  sufficient  range  of  saturation  tem- 
peratures and  not  exceed  the  35  psia  limit  on  the  heat  pipe, 
distilled  water  ice  was  placed  in  the  evaporator  and  a  low 
vacuum  was  drawn.   The  filling  procedure  was  as  follows: 

a.  System  was  cooled  by  circulating  cooling  water 
through  the  condenser  nozzles  and  by  inserting  distilled 
water  ice  in  the  evaporator. 

b.  Six  distilled  water  ice  cubes  (approximately 
150  ml)  were  subcooled  by  flashing  them  with  liquid  N2  and 
placed  into  the  evaporator  with  tongs. 
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c.  Glass  windows  were  put  into  place  and  the  flange 
was  tightened  with  a  30  in-lb.  torque. 

d.  The  vacuum  was  drawn  through  the  glass  tube 
connected  to  a  mechanical  vacuum  pump. 

e.  While  the  vacuum  was  being  pumped,  the  glass 
tube  was  sealed  using  a  gas  flame. 

2 .   Operational  Procedure 

Having  completed  the  filling  procedure,  the  heat 
pipe  is  ready  to  run  at  a  pre-determined  PPM.   About  10 
drops  of  lubricating  oil  are  added  to  the  oil  fill  tube 
at  the  top  of  the  large  pillow  block.   The  cooling  water 
lines  are  turned  on  so  that  the  rotameter  is  at  approxi- 
mately 20%  flow. 

The  drive  motor  is  turned  on  and  the  heat  pipe  is 
slowly  brought  up  to  the  desired  RPM.   Once  a  liquid 
annulus  has  formed  in  the  evaporator  and  the  desired  RPM 
has  been  reached,  the  DC  motor  generator  is  turned  on  and 
power  is  added  to  the  evaporator.   Having  reached  thermal 
equilibrium  for  a  particular  power  setting,  data  was  taken 
by  the  instrumentation  described  earlier. 

Shut-down  is  accomplished  by  first  shutting  off  the 
current  to  the  heater,  and  then  decreasing  the  RPM  by  use 
of  the  remote  control  unit. 

C.   EXPERIMENTAL  RESULTS 

The  heat  pipe  was  run  at  702  RPM  and  1404  RPM.   The  run 
of  702  RPM  was  chosen  because  it  was  the  lowest  RPM  at  which 
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centrifugal  forces  dominated  the  gravitational  effects. 
Power  inputs  to  the  heat  pipe  ranged  from  1  kW  to  9  kW. 

1 .   Preliminary  Results 

Using  distilled  water  as  the  working  fluid,  the  heat 
pipe  was  rotated  at  70  2  RPM  to  obtain  the  heat  flow  through 
the  evaporator  wall  for  a  particular  power  setting.   The 
method  used  to  compute  0     ,  the  heat  flow  through  the 
evaporator  wall,  and  the  power  supplied  to  the  heater  coils 
is  shown  in  Appendix  B.   Table  II  shows  a  comparison  between 
the  electrical  power  into  the  heater  coils  and  the  heat  flow 
through  the  evaporator  wall  for  thcit  particular  power 
setting. 

The  data  taken  from  the  four  thermocouples  in  the 
evaporator  wall  is  obviously  incorrect  since  part  of  the 
data  shows  that  the  heat  flow  through  the  evaporator  wall 
is  higher  than  the  electrical  heat  in  for  that  particular 
power  setting.   In  other  cases  it  appears  to  be  low. 

There  are  two  explanations  for  the  incorrect  measure- 
ment of  the  heat  flow  through  the  evaporator  wall.   The 
first  is  the  fact  that  the  reading  of  thermocouple  #2 
appeared  erratic.   When  the  heat  pipe  is  not  rotating, 
thermocouple  #2  appears  to  be  reading  correctly,  but  v/hen 
the  heat  pipe  rotates  it  begins  to  read  inadequately.   It 
was  concluded  that  the  problem  existed  in  the  slip  ring  and 
may  be  due  to  insufficient  mercury  in  the  slip  ring.   The 
second  explanation  pertains  to  the  reason  why  the  heat  flow 
through  the  evaporator  wall  is  higher  than  the  heat  into 
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TABLE  II 
POWER  IN  AND  HEAT  FLOW  THROUGH  EVAPORATOR  WALL 


POWER  IN 


HEAT  FLOW  IN 
EVAPORATOR  WALL 
(BTU/HR)  (BTU/HP) 


3,345  4,059 

6,712  10,315 

9,963  12,848 

13,325  7,610 

16,988  9,126 

20,751  24,041 

24,031  13,437 

27,278  31,986 


the  heater  coils  for  a  particular  power  setting.   It  was 
concluded  that  the  outside  edge  coils  on  each  side  of  the 
heating  surface  were  conducting  axially.   As  a  consequence, 
there  wasn't  a  uniform  temperature  distribution  axially 
along  the  heating  surface  as  assumed.   Therefore,  the  actual 
heating  surface  length  was  less  than  the  assumed  length  in 
the  calculation  of  the  heat  flow  and  would  account  for  the 
high  values  of  the  heat  flow  obtained  in  Table  II. 

At  this  point,  it  was  decided  to  obtain  the  heat  flow 
out  of  the  heat  pipe  by  performing  an  energy  balance  on  the 
condenser  cooling  section.   The  coolant  flow  rate  was 
determined  by  a  flow  rotameter,  and  its  mean  temperature  was 
measured  before  and  after  spraying,  by  quartz  thermometers 
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mounted  in  mixing  tubes  located  above  and  below  the  conden- 
ser box.   Having  this  information,  along  with  the  coolant's 
specific  heat,  the  heat  flow  was  determined  by  a  simple 
energy  balance.   In  order  to  make  this  method  more  accurate , 
the  inside  of  the  condenser  box  was  insulated  with  silicone 
foam  rubber  and  the  plastic  coolant  lines  were  insulated 
with  fiberglass.   Appendix  B  contains  a  table  showing  the 
heat  flow  out  of  the  condenser  at  different  values  of  the 
heat  pipe  saturation  temperature. 

2 .   Visual  Observations 

Using  a  strobe  light  visual  observations  were  made 
through  the  view  ports.   At  702  RPM,  the  light  provided  by 
the  strobe  light  was  insufficient.   To  increase  the  light, 
the  strobe  light's  frequency  was  increased  to  three  times 
the  RPM  of  the  heat  pipe.   As  a  consequence,  the  observations 
made  at  702  RPM  were  not  clear.   The  observation  of  the 
liquid  vapor  interface  in  the  evaporator  at  low  saturation 
pressures  indicated  violent  boiling  with  large  vapor 
bubbles  breaking  through ' the  liquid  vapor  interface  at 
702  RPM.   At  higher  saturation  pressures,  the  boiling  was 
more  controlled  with  small  bubbles  constrained  to  the 
surface  of  the  liquid.   At  the  same  saturation  pressure,  the 
boiling  was  less  violent  at  1404  RPM. 

The  visual  observation  in  the  condenser  indicated  that 
there  were  ripples  and  random  break-up  of  the  condensate 
film,  exposing  parts  of  the  condenser  surface  to  the  vapor. 
Furthermore,  the  condensate  swirled  and  washed  the  condenser 
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wall  at  a  higher  RPM  than  the  wall  itself.   It  is  doubtful 
that  the  above  observed  effects  were  due  to  condensate 
turbulence.   According  to  Rohsenow  and  Choi  [3],  in  the 
presence  of  interfacial  shear,  turbulent  condensate  flow 
occurs  for  a  Reynolds  number  between  100  and  1800.   For  the 
low  vapor  velocities  which  occurred  during  the  observed 
experimental  runs ,  the  calculated  transition  Reynolds  number 
in  the  condensate  was  approximately  100  0.   However,  the 
highest  Re  calculated  for  the  condensate  was  only  291.   It 
was  concluded,  therefore,  that  the  swirling  and  random  film 
break-up  were  not  caused  by  condensate  turbulence  but  were 
probably  caused  by  vapor  interaction  with  the  condensate. 

Except  for  the  1  kW  case,  the  Reynolds  number  calculated 
for  the  vapor  flow  for  all  of  the  power  inputs  was  between 
2000  and  9000  (for  1  kW,  Re  =  972) .   Since  flow  in  a  smooth 
tube  is  turbulent  for  Reynolds  numbers  between  2000  and 
50,000  [3],  it  was  concluded  that  the  vapor  flow  was  tur- 
bulent except  for  the  1  kW  case.   The  swirling  of  the 
condensate  was  probably  caused  by  the  higher  velocity  vapor 
striking  the  condensate. 
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III.   DISCUSSION  OF  RFSULTS 

The  graphs  in  Figures  9  and  10  show  the  experimental 
and  theoretical  heat  transfer  rates  as  a  function  of  the 
saturation  temperatures  (corresponding  to  an  internal  vapor 
pressure) . 

In  order  to  find  the  theoretical  values  of  the  heat 

transfer  rate,  the  value  h,  the  heat  transfer  coefficient 

was  needed.   The  method  used  in  determining  h  is  shown  in 

Appendix  B.   Because  of  experimental  uncertainties,  the 

values  of  h  .   =  1000  and  h     =  1800  were  used  since  these 
mm  max 

values  were  very  close  to  the  average  of  the  minimum  and 
maximum  experimental  values  of  h  determined  for  the  dif- 
ferent operating  conditions. 

As  observed  in  Figure  9 ,  the  theoretical  values  of  Q  . 

and  Q     are  less  than  the  experimental  value  of  Q  at 
max 

702  RPM.   In  the  development  of  the  analytical  solution  it 
was  assumed  that  smooth  fi'lm  condensation  occurred  in  the 
condenser.   The  visual  observations  described  in  the  experi- 
mental results  indicated  that  a  smooth  condensate  film  did 
not  exist  for  the  operating  conditions  used  in  these 
experiments  but  rather  that  the  condensate  film  swirled  and 
washed  the  condenser  wall  and  appeared  to  break  up  in  a 
random  manner  exposing  part  of  the  condenser  surface  to  the 
vapor.   The  swirling  effect  and  the  resulting  increased 
convection  in  the  condensate  increased  the  heat  transfer 
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rate.   It  was  also  possible  that  the  ripples  and  film 
break-up  decreased  the  overall  film  thickness,  reducing  the 
thermal  resistance,  and  thereby  increasing  the  heat  transfer 
rate. 

Figure  10  shows  that  the  analytical  and  experimental 
results  compare  more  favorably  for  the  1404  RPM  run.   There 
are  two  probable  explanations.   First  of  all,  the  higher 
RPM  could  have  flattened  the  film  causing  it  to  approach  a 
smoother  film  than  in  the  702  PPM  run.   If  this  were  the 
case,  the  smooth  film  condensation  assumptions  would  be 
better,  and  the  analytical  solution  would  be  closer  to  the 
experimental  results.   Although  the  visual  observations 
seemed  to  indicate  that  there  were  no  non-condensibles 
present,  the  presence  of  non-condensibles  in  the  condenser 
could  be  the  second  explanation.   If  non-condensibles  were 
present,  the  experimental  heat  transfer  curve  would  be 
shifted  to  the  right  since  the  heat  transfer  rate  would 
require  higher  vapor  pressures  (and,  therefore,  higher 
saturation  temperatures). 

If  the  inlet  coolant  temperature  were  held  constant, 
the  heat  transfer  rate  should  increase  with  increased  RPM 
for  the  same  saturation  temperature.   For  example,  the 

analytical  solution  for  Q  .   at  1404  RPM  at  apnroximately 

"mm  ■ 

the  same  T^  used  at  702  RPM  resulted  in  higher  heat  transfer 
rates  as  shown  by  the  dashed  plot  in  Figure  10.   Note  that 
the  experimental  data  at  1404  RPM  does  not  appear  to  show 
an  increase  in  heat  transfer  rates  from  the  70  2  RPM  run. 
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Figure  10.   Heat  transfer  rate  vs  saturation  temperature. 
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This  is  due  to  the  fact  that  the  1404  RPM  run  was  conducted 
at  a  higher  T^  (7°  difference)  than  the  702  PPM  run.   This 
increase  in  Tm  increased  the  saturation  temperature  by 
approximately  the  same  amount  such  that  the  experimental 
heat  transfer  rate  data  points  were  shifted  to  the  right. 
At  the  same  saturation  temperature,  therefore,  the  heat 
transfer  rate  was  lower  for  the  higher  Tw  even  though  it 
was  at  higher  RPM. 


49 


IV.   CONCLUSIONS  AND  RFCOMMFNDATIONS 

A.   CONCLUSIONS 

1.  The  rotating  non-capillary  heat  pipe  is  an  effective 
heat  transfer  device. 

2.  Smooth  film  condensation  was  not  observed  during  the 
experimental  runs. 

3.  The  approximate  numerical  solution  conservatively 
predicts  the  heat  transfer  rate  with  a  deviation  of  18%  at 
702  RPM.   The  deviation  of  the  numerical  solution  decreases 
to  5%  for  1404  RPM. 

4.  The  performance  of  the  heat  pipe  is  strongly 
dependent  on  the  method  used  to  cool  the  condenser. 

B.   RECOMMENDATIONS 

1.  Run  the  heat  pipe  at  higher  RPM  to  determine 
whether  increased  RPM  provides  better  agreement  between  the 
analytical  and  experimental  results. 

2.  Increase  the  instrumentation  of  the  condenser  to 
give  more  accurate  data  for  the  experimental  heat  transfer 
rate  out  of  the  condenser. 

3.  Revise  the  analytical  solution  to  include  the  ripples 
and  film  break-up  in  the  condenser. 

4.  Obtain  experimental  results  for  the  effect  of  non- 
condensibles  in  the  heat  pipe. 

5.  Use  different  methods  of  cooling  the  condenser  i 
by  using  different  spray  nozzles. 
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6.  Use  a  condenser  section  with  a  half -cone  angle 
greater  than  1°  to  obtain  experimental  data  for  the  effect 
of  larger  cone  angles  on  the  heat  transfer  rate. 

7.  Use  working  fluids  with  boiling  points  lower  than 
water  in  order  to  find  optimum  heat  transfer  rates  at  lower 
saturation  temperatures. 
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APPENDIX  A 


CALIBRATION 


1.  Calibration  of  Pressure  Transducer.   The  pressure  trans- 
ducer consists  of  a  temperature-compensated,  semiconductor 
strain  gage  Wheatstone  bridge  mounted  on  a  1/4-inch  diameter 
diaphragm.   It  was  calibrated  against  a  Wallace-Tiernan 
standard  pressure  gauge  (with  a  known  accuracy  of  0.03  ± 
0.01  psia)  for  pressures  greater  than  atmospheric,  and 
against  a  mercury  manometer  (with  an  uncertainty  of  ±.02 
psia)  for  pressures  less  than  atmospheric.   The  transducer 
was  calibrated  at  different  operating  temperatures  by 
placing  it  within  a  controlled  temperature  oven  with  a 
quartz  crystal  thermometer  (with  an  uncertainty  better  than 
±0.1°F).   Reproducible  data  was  obtained  from  0-30  psia 

and  from  65-222°F,  and  the  transducer  sensitivity  was 
measured  to  be  3.54  ±0.06  mv/psia.   Although  the  transducer 
Wheatstone  bridge  was  temperature  compensated,  a  non-linear 
zero  shift  with  temperature  was  measured  and  incorporated 
into  a  series  of  calibration  curves. 

2.  Calibration  of  Thermocouples.   Nine  30  gauge,  Kapton 
insulated,  copper-constantan  thermocouples  were  carefully 
manufactured  at  the  Naval  Postgraduate  School.   They  were 
calibrated  up  to  384°F  using  the  same  temperature  controlled 
oven  and  quartz  crystal  thermometer  mentioned  above.   In 
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addition,  the  following  fixed  point  calibrations  were  made: 

a.  Triple  point  of  water  (32.0°F) 

b.  Boiling  point  of  water  (212. 0°F  at  14.696  psia) 

c.  Melting  point  of  tin  (449. 4°F) 

d.  Melting  point  of  lead  (621. 3°F) 

There  were  a  total  of  nine  calibration  data  points 
taken  over  the  range  of  32-621.3°F.   Using  the  NBS  table 
of  calibration  data  for  copper-constantan  thermocouples  as 
a  zero  reference,  a  deviation  curve  of  the  error  (mv)  vs. 
the  thermocouple  reading  (mv)  was  plotted.   This  curve  (with 
an  uncertainty  of   0.1%)  was  then  used  as  a  correction 
factor  for  each  data  point  taken  during  the  experiment 
before  entering  the  NBS  table  for  the  temperature.   Recali- 
bration  data  was  taken  at  a  later  time  and  the  repeatability 
of  all  nine  thermocouples  was  excellent.   Using  the  boiling 
point  of  water  and  the  melting  point  of  tin,  the  deviation 
was  within  0.00  5  mv  of  the  previous  data  points. 

3.   Calibration  of  Rotameter.   The  method  used  in  the  cali- 
bration of  the  flow  rotameter  consisted  of  using  a  Toledo 
scale  to  measure  the  mass  of  the  coolant  flowing  into  a 
bucket  placed  on  the  scale.   As  coolant  flowed  into  the 
bucket,  an  arbitrary  starting  point  on  the  scale  was  taken 
and  a  stop  watch  was  used  to  measure  the  time  it  took  for 
the  weight  to  reach  20  pounds.   This  procedure  was  repeated 
several  times  for  each  per  cent  flow  rate  used  on  the 
rotameter.   Data  points  were  obtained  for  14,  16,  18  and 
19.7  per  cent  flow  rates.   From  this  data  a  linear  calibra- 
tion curve  was  obtained. 
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DATA  RFDUCTION 


1.   Calculation  of  Q 


EVAP 


The  heat  flow  through  the  evaporator  wall  may  be  reduced 
to  a  one  dimensional  steady  state  problem  if  it  is  assumed 
that  there  is  no  axial  heat  conduction  and  steady  state 
conditions  exist.   Under  these  conditions,  the  well  known 
conduction  equation  in  cylindrical  coordinates  may  be 
written  [1] : 


1  _d_ 
r  dr 


k(T)  r  ^ 

Ui. 


0 


k(T)  r  |^-  =  B  =  Constant 


(a) 


where : 


k(T) 


k   +  m(T  -  T  ) 
o         o 


(b) 


A  plot  of  k  vs .  T  was  made  based  on  data  taken  from  the 
Metals  Handbook  [5]  for  the  thermal  conductivity  at  different 
temperatures.   From  this  plot  the  slope,  m,  and  k   were 
determined . 

Substituting  equation  (b)  into  equation  (a)  and  inte- 
grating yields : 


F(T) 


k   +  m 
o 


T  =  B  In  r  +  c. 
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Using  a  sub-routine  for  the  method  of  least  squares 
the  IBM  360  computer  was  used  to  find  the  slope,  B,  of 
F(T)  vs.  In  r.   From  Fourier's  law: 


QEVAP  =  L  2 


*[*Ct,  r«] 


or 


where 


QEVAP  "  L  2*  B 


L  =  the  effective  heating  surface  length  in  the 
evaporator 


EXAMPLE:   T     =  153. 7°F,  702  RPM 


QEVAP  =  L  2tt  B 


L  =  .328  (ft. ) 

B  =  4430.  (BTU/hr-  ft.) 

QEVAP  =  9,126  (BTU/hr) 
2.   Calculation  of  POWER  IN 


PT.7  =  Voltage  x  Current 
IN 


3.   Calculation  of  Experimental  QnTTm 


where 


QOUT  "  m  Cp(Tl"T2) 


m  =  mass  flow  rate  (lbm/hr) 


C   =  specific  heat  (BTU/lbm  °F) 
T,  =  inlet  coolant  temperature  (°F) 
T„  =  exit  coolant  temperature  (°F) 
Q0UT  =  heat  flow  (BTU/hr) 


EXAMPLE:   TSAT  =  153. 72°F,  702  RPM 


m  =  797  ±  l: 


Cp  =  1.0 

(Tx  -  T2)  =  16.04  ±0.31 
Q0UT  =  12,784  1  541 


Calculation  of  h, outside  condenser  wall  heat  transfer 
coefficient 


°OUT=    hVTw-T~> 


h    =    0oOT/Ac(Tw-,T-) 


where 


2 
A   =  cooling  surface  area  (ft  ) 

T  =  outside  wall  temperature  of  condenser  (°F) 

Tot  =  coolant  inlet  temperature  (°F) 

h   =  heat  transfer  coefficient  (BTU/hr-f t2-°F) 

EXAMPLE:   TSAT  =  153.72°,  702  RPM 

QQUT  =  12,784  1  541 

A   =  .34 
c 

(T   -T  )  =  28.44  ±  6.7 
w     °° 

h  =  1405  ±  386 
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5.   Tables  of  Reduced  Data 


T,  -  T 

T   -  T 

Q 

h 

TSAT 

m 

1     2 

W     °° 

^OUT 

AVG  7 
(BTU/HR-FT  -°F) 

(°F) 

(LBM/HR) 

(°F) 

(°F) 

(BTU/HR) 

70  2  RPM 

93.4 

835±23 

3.061.11 

6.611  1.8 

2,5551165 

12421411 

107.3 

832±23 

5.921.25 

11.1  i  2.6 

4,9251350 

14031423 

122.8 

804±20 

9.231.31 

17.2  l  2.7 

7,4211440 

13051280 

137.8 

823±24 

12.5  1.54 

22.9  l  4.6 

10,2961757 

13881375 

153.8 

797±18 

16.0  1.31 

28.4  l  6.7 

12,7841541 

14051386 

172.1 

764±19 

20.5  1.25 

34.1  i  7.2 

15,6471585 

14171350 

188.3 

764±19 

23.6  1.31 

38.9  H0.3 

17,9921690 

15151465 

204.4 

863119 

22.6  l .34 

37.9  l  8.5 

19,5381730 

16491438 

14  0  4  RPM 

98.4 

870123 

3.021.04 

1.0   11.1 

2,6271102 

12251408 

115.3 

850124 

5.961.18 

12.6   15.2 

5,0661300 

14451664 

141.8 

825122 

12.2  1.21 

21.7   t4  .  0 

10  ,0  821438 

1 4 .1 7 1  3  2  C 

174.3 

790120 

18.9  1.53 

31.7   16.5 

14.9391783 

1456±371 

205.6 

786127 

25.4  1.41 

36.1   13.9 

19,96411011 

16431259 

224.3 

786127 

29.4  1.68 

45.6   19.0 

23,13211347 

15631395 
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